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Abstract 

We present the Kormendy and mass-size relations for early-type galaxies (ETGs) as a function of 
environment at z ~ 1.3. Our sample includes 76 visually classified ETGs with masses 10^° < M/Mq < 
10"•^ selected in the Lynx supercluster and in the GOODS/CDF-S field, 31 ETGs in clusters, 18 in 
groups and 27 in the field, all with multi- wavelength photometry and HST/ ACS observations. The 
Kormendy relation, in place at z ~ 1.3, does not depend on the environment. The mass-size relation 
reveals that ETGs overall appear to be more compact in denser environments: cluster ETGs have 
sizes on average around 30-50% smaller than those of the local universe, and a distribution with a 
smaller scatter, whereas field ETGs show a mass-size relation with a similar distribution than the 
local one. Our results imply that (1) the mass-size relation in the field did not evolve overall from 
z ~ 1.3 to present; this is interesting and in contrast to the trend found at higher masses from previous 
works; (2) in denser environments, cither ETGs have increased their size by 30-50%, on average, and 
spread their distributions, or more ETGs have been formed within the dense environment from not 
ETG progenitors or larger galaxies have been accreted to a pristine compact population to reproduce 
the mass-size relation observed in the local Universe. Our results are driven by galaxies with masses 
M < 2 X 10"Mo and those with masses M ^ 10^"'^ Mq follow the same trends that the entire sample. 
Following Valentinuzzi et al. definition of superdense ETGs, around 35-45% of our cluster sample is 
made of superdense ETGs. 

Subject headings: galaxies: clusters: individual (RX J0849-I-4452, RX J0848-I-4453) - galaxies: ellip- 
tical and lenticular - galaxies: evolution - galaxies: formation - galaxies: high- 
redshift - galaxies: fundamental parameters 
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1. INTRODUCTION 

In recent years, studies have unveiled the existence at 
z ~ 1-2 of a population of massive sphero idal galaxies 
with sm all size, hence c alled c on ipact fe. g..lDaddi et al " 
^005t iTruiillo et all [20061 [20071: iBuitrago et al 



2008t iCimatti et all l2008t Ivan der Wei et a.1 1 



van Dokkum et al. 120081 Darnianov et alT 



baracco ct al.' '2009; 'Newman ct al.' 2010'; 'Rc ttura et al.l 
2010,;.Saracco ct a l. 2010; Strazzullo ct al. .2010i and also 



references therein) . When comparing those high redshift 
galaxies with local ones of similar mass, it appears that 
th eir sizes are smaller by a factor of ^2-3 and up to 
5 (|van Dokkum et all |2008| ) . The general view is that 
the compactness increa ses with redshift, mass and the 
level of quiescence (e. g.. ITruiillo et al.l[2007l : iFranx et all 
120081 : IWilliams et al.ll2010[ )._ Despite potential selection 
biases affecting the comparison of high vs low redshift 
samples (see hereafter), which might affect conclusions 
on the evolution in size, the existence of a significant 
number of compact galaxies at high redshift is firmly 
established. [3- 

The presence of compact ETGs in the local Universe is 
still debatable. Apparent disagreements may come from 
the different definitions for a compact galaxy (i.e. the dif- 
ferent mass and size criteria chosen to define a galaxy as 
compact). For example, on the o ne hand, the analy sis of 
Sloan Digital Sky Survey (SDSS; lYork et al.l[200(if) sam- 
ples r eveals that a neglig ible fraction of galaxies are com- 
pact (ITruiillo et al.|[2009[ ). even when taking into account 
the possible incompleteness due to the SDSS spe ctro- 
scopic target se lection algorithm (iTavlor et al.ll2010D . On 
the other hand, IValentinuzzi et alT i 2010af ) studied ETGs 
in local clusters and found that a significant fraction of 
their sample is made of comp act objects. When com- 
pared to high redshift samples (jSaracco et al.|[2009[ ). the 
number density of compact galaxies at 1 < z < 2 is con- 
sistent with that found in this last work and consistent 
with a lack of evolution in size (see also Shankar ct all 
[20T0I : IBernardi et al.llMol) . 

The formation of compact galaxies might be a 
consequence of mergers of gas- r ich subunits at high 
redshift Ce.g.. IKhochfar fc Silkl l2006t iHopkins et ahl 
2009bt iWuvts et al.l I2010D and /or cold flows (e.g. 



Bournaud et al.l 120111 ). resulting in an intense star- 
bust and compact quiescent remnant due to highly 
dissipative processes. This is in agreement with 
observations showing that the gas fra ction of star- 
formi ng galaxies increases with redshift (jHopkins et al.l 
|2010[ ). Sub- millimeter galaxies have been suggested as 
Tomising cand idates for compact galax ies precursors 
Granato et £lll2006l : ICi^atti et al.l [20081 ). The picture 
concerning the subsequent evolution of compact galax- 
ies down to z = is more difficult to draw. The com- 
parison of high redshift to local samples may be af- 
fected by two selection biases: age selection bias against 
youiig galaxies in high redshift samples fe.g. lSaglia et al.l 
I2OIOI : IValentinuzzi et al.ll2010al ) a nd progenitor bias due 
to ni o rphological evolution (e.g.. Ivan Dokkum fc Frarp3 
[2OOI iKavirai et all [2001 IValentinuzzi et al.l l2010bi r 
Within this context, it is still unclear which part of 



(e.[ 

|oi 



We underline that n ot all high redshift ETGs are compact 



McGrath et a l.'2008f:'Saracco ct al.l200gi : IMancini et al.l2010l : 
et al.ll20l0., : .Saracco ct al. 20lJ) 



the galaxy population went through evolution and which 
mechanism contributed to it. If the compact galaxy 
population requires evo lution, one efficient process may 
be m inor dry mergers (jNaab et al.ll2009t iShankar et al.l 
l2011h : through the accretion of gas-poor satellites, a 
compact galaxy will increase significantly its size with 
a limited increase of its mass and no star formation. In 
this scenario, the accreted material will extend the outer 
parts of the compact galaxy, leaving its core unchanged. 
This is in remarkable agreement with observations: local 
elliptical galaxies have in their core regions surface stel- 
lar density profiles similar to those of high-redshift com- 
pact galaxies (e.g.. iBezanson et al. [ |2009l : [Hopkins et al.l 
I2009at [ van Dokkum et al.ll2010[ ). Another proposed sce- 
nario for size evolution of compact spheroids is expansion 
consequent to substantial mass losses due to, e. g., stellar 
winds and/or quasar feedback (|Fan et al.ll2008l ). 

To go deeper in understanding these mechanisms, it 
is useful to study the mass-size relation as a function of 
environment. Until now, few studies have covered the 
full range of environment when studying the mass-size 
relation at z > 1. In the local universe, iMaltbv et al.l 
(|2010( ) have found that the mass-size relation does not 
depend on the environment for ETGs. Most of the cur- 
rent z > 1 studies tho ugh r ely on field sampl e s, exc ept for 
iRettura et all (I2010D and [Strazzullo et al l (I2010D. who 
studied clusters at z ^ 1.2-1.4. Onlv IRettura et al.l 
(|2010f ) compared field and cluster ETGs at z ~ 1.2 and 
find that galaxies from different environments lie on the 
sam e relations. 

In [Raichoor et al.l (|2011l Rll herafter), we presented 
a unique homogeneous sample of ETGs probing cluster, 
group and field environments at z 1.3. Our study relies 
on high-qualit y multi- wa velength data covering the Lynx 
supercluster jStanford ct al. 1997; Rosati ct al.. 19991: 
Nakata et aTll2005HMei'erar]l2006bLl201lURettura et al.l 
20111) . a structure at z = 1.26 made of two clusters 
and at least three groups. From our spectroscopic runs 
on the groups, we obtained average spectroscopic red- 
shift z — 1.262 ± 0.007 (Group 1; from 9 members), 
z = 1.260 ± 0.006 (Group 2; from 7 memb ers), z = 
1.263 ± 0.005 (Group 3; from 9 members) (|Mei et al.l 
[20 111 ). Group X-ray emission g ives masses less or around 
5 X IO^^Mq (|Mei et al.l[2011l ). Group 2 and 3 appear 
to be spatially separated (as from our Fricnd-of-Friend 
algorithm) from the two clusters, while Group 1 is spa- 
tially connected to the Lynx W cluster. We consider 
it as a separate group though, because its center is at 
1.1 X r2oo from the center of the cluster and it extends 
to 2 X r2oo, with an area of very low density between 
0.5-1x^200- It might be close to merging to Lynx W, 
or in the merger pr ocess. For further details please refer 
to lMei et al.l (|201ll ). Our groups belong all to the Lynx 
supercluster, and are not isolated. They do not show 
peculiar densities, or masses to differentiate them from 
isolated groups. A more extended analysis of superclus- 
ter groups as compared to isolated groups at the same 
redshift would help us understanding if the properties of 
their galaxies might be different. At the moment, we do 
not have elements to suggest it. 

In this paper, we use the Rll sample to study the in- 
fluence of environment on the structural parameters of 
ETGs at high redshift as a function of mass and envi- 
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ronmcnt. The estimates of ETG sizes from HST/ ACS 
images combined with the photometry and the stellar 
population parameters determined in Rll allow us to 
build the two key relations to study structural parame- 
ters of ETGs: the IKormendvl (|1977D relation (KR) and 
the mass-size relation (MSR). 

The plan of this paper is as follows. In §2, we present 
the observations, the sample selection and the SED fit- 
ting method used to estimate ages and masses. In §3, we 
describe our estimation of the ETG structural parame- 
ters. In §4, we study the Kormcndy relation and in §5 
the mass-size relation. We then present our conclusions 
in §6. 

We adopt a standard cosmology with Hq = 70 km s~^ 
Mpc-\ n^n = 0.30 and VIa = 0.70. All magnitudes are 
in the AB system. Unless otherwise stated, all stellar 
masses are computed with a lSalpeteil (|1955D Initial Mass 
Function (IMF). We choose as our rest-frame reference 
the Coma cluster (zq = 0.023). 

2. OBSERVATIONS, SAMPLE SELECTION. PHOTOMETRY 
AND SED FITTING 

This work relies on optical and infrared (0.6-4.5 
/xm) images of the Lynx supercluster and of the 
Great Observatories Origins Deep Survey (GOODS; 
iGiavalisco et al.l |2004[ ) observa tions of the Chandra 
Deep Field South CCDF-S: IGiavalisco et all [200l 
INonino et sl\ I2009t IRetzlaff et al.l 120101: Dickinson et 
al., in preparation). The observations, the sam- 
ple selection, the photometry and the age and stel- 
lar mass estimation are presented in Rll and we 
briefly summarize them here; please refer to Rll for 
more details. The images cover seven bandpasses: R 
(Keck/LRIS for the Lynx clusters, Palomar/COSMIC 
for the Lynx groups, VLT/VIMOS for the CDF-S), 
HST /ACS F775W and F850LP - hereafter ^775 and zgso, 
J/K, (KPNO/FLAMINGOS for the Lynx clusters and 
groups, VLT/ISAAC for the CDF-S), Spitzer /IRAC chl 
and ch2 - hereafter [3.6/im] and [4.5/im]. The sample 
of Rll consists of 79 ETGs (31 in the Lynx clusters, 21 
in the Lynx groups and 27 in the CDF-S) selected in 
redshift (0.92 < Zphot < 1-36 for the Lynx ETGs and 
1.1 < Zspec < 1-4 with (zspec) = 1-239 ± 0.082 for the 
CDF-S), in magnitude (21 < zgso (AB) < 24) and in mor- 
phology (E/SO types based on visual ins pection of zh^o- 
band of H ST /ACS images according to iPostman et al.l 
(12005) and iMei et all (|2011l ) classification). We verified 
that the magnitude cut zg^Q > 21 does not exclude any 
galaxy satisfying the Zphot or z^pec selection criteria; thus 
we can relax the magnitude cut to 2350 (AB) < 24 with- 
out affecting the sample. ETGs belo nging to th e Lynx 
clusters and groups are identified in iMei et al.l (1201 ID 
by a Friend-O f- Friend algorithm (F OF. lGeller fc Huchral 
119831 see also IPostman et all 120051 ) with a linking scale 
corresponding to a local distance of 0.54 Mpc, normalized 
to z = 1.26 and to ou r magnitude range (jPostman et al.l 
[20051 IMei et"aI1[20Tlh . We also verified that the selected 
CDF-S ETGs are field ETGs, i.e. that they do not be- 
long to already identified structures (see Rll). 

At zgso = 24 mag. Lynx samples are complete and our 
CDF-S sample is more than 70% complete (see Rll). 
The Lynx cluster, group, and CDF-S field samples have 
similar spectral coverage and are almost complete at 
•^850 = 24 mag, thus providing a homogeneous and consis- 



tent sample. Since the publication of Rll, spectroscopic 
observations revealed that three ETGs from our Group 
2 sample were outhers (ID = 939, 1791, 2519). We thus 
remove those three ETGs from our sample, obtaining a 
final sample of 76 ETGs (31 in the Lynx clusters, 18 
in the Lynx groups and 27 in the CDF-S). The removal 
of those three outliers does not affect significantly any 
of the results presented in Rll. Our sample has spec- 
troscopic redshifts for 20/31 ETGs in the clusters, 8/18 
ETGs in the groups (|Mei et al.ll20TT[ ) and 27/27 ETGs 
in the field. 

We performed photometry in circular apertures with 
1.5" radius and derive a multi- wavelength photomet- 
ric catalog with total magnitudes, determined using PSF 
growth curves. We estimated stellar masses and stel- 
lar population ages by fitti ng the SED with d ifferen t 
stellar po p ulatio n models (jBruzual fc Charlotl (|2003| ). 
MarastonI (I2005D, and a n updated version [CB07] of 



Bruzual fc Charlotl (|2003[ ) that implements a new mod- 
eling of the TP-AGB phase). We hereafter refer to 
those models as BC03, M A05 and CB07, respectively. 
For SED fitting we used a iSalpeted (|1955f l IMF, solar 
metallieity, exponentially declining star-formation histo- 
ries ip{t) oc e~*/'^ with a characteristic time 0.1 < SFH 
T (Gyr) < 5, and no dust. Our stellar mass is the mass 
locked into stars, including stellar remnant^ and our 
age is star-formation weighted age. A detailed discus- 
sion of different choices of parameters can be found in 
Rll. 

3. SIZE ESTIMATION 
3.1. Method 

In this section, we describe our methodology to de- 
rive the size of our ETGs. Morphological parameters 
are usually estimated in the rest-frame i?-band: we de- 
rive them from the HST/ ACS zgso band image, the clos- 
est to the rest-frame S-band in our sample. To fit 
the observed two-dimensional surface brightness distri- 
butions to a mod el, we use the software Galfit (v3.0.2, 
iPeng et "all l2002f). which has be en shown to giv e reli- 
able results (jHaussler et al.l 120071 ). We assumed a lSersid 
(p68h 



profile: 



I{r) = /e X exp{-bn[{r/re) 



l/n _ 



1]}, 



(1) 



where /(r) is the surface brightness at r, Ig is the surface 
brightness at the effective radius r^, which is the radius 
which encloses half of the emitted light. In the fit, Gal- 
fit convolves the model with a provided PSF: our PSF 
stamp is built from real isolated unsaturated stars, by 
first normalizing them and then taking the median value 
for each pixel (the same as the one used in Rll, see this 
paper for more details). It has been shown in the liter- 
ature that considering different stars for the PSF leads 
to mi nor changes in the size estimate (e.g. iTrujillo et al.l 
|2007| ). Galfit outputs the semi-major axis of the 
projected elliptical isophotc containing half of the total 
light and the axis ratio b/a. Throughout this work, we 
use Re to denote the circularized effective radius defined 

Using the nomenclature given by the authors; column 7 of 
*.4color files for BC03/CB07 models and "MA* total" for MA05 
models 
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by: 

Re = re X v/Va- (2) 

For each object, we create a square stamp from the 
ACS image centered on the galaxy. According to our 
tests (using stamp size of 2.5 x ri, 5 x ri and 10 x ri), 
the fit is stable for a stamp size of 5 x ri , where ri de- 
note s thelKronI (119801) radius , as determined by SExtrac- 
tor (jBertin &: Arnoutill996[) . We simultaneously fit the 
selected ETG along with any objects closer than 2.5"and 
use SExtractor segmentation maps to mask the other ob- 
jects. During the fit, we let as free parameters the posi- 
tion {x,y), the total magnitude zgso, the effective radius 
Te, the axis ratio 6/a, the Sersic index n and the posi- 
tion angle pa. While we use SExtractor outputs as ini- 
tial guess for {x,y), zgso, r^, b/a and pa, we set the initial 
Sersic index n to 2.5. As advised in the Galfit homepage, 
no boundary constraint on the Sersic index is provided 
during the fit, so that the minimization algorithm can 
run properly. In order to reduce the number of free pa- 
rameters and improve the quality of the fit, we fix the 
sky value. For sky estimation, we create a larger stamp 
(20"x20") centered on the ETG, we mask the objects 
with ellipses (taking SExtractor's outputs and increasing 
the linear size by a factor 5) and take the median value 
of the remaining pixels. This conservative approach for 
masking objects ensures that there is negligible residual 
light from the objects in the sky area, while keeping a 
large enough number of pixels. 

For five ETGs of our sample 6%), our fits do not 
provide satisfactory results: either the output parame- 
ters are unphysical {n > 10, small Re)-, or the residuals 
are unsatisfactory (two ETGs). For those five ETGs, 
we consider the structural parameter estimates as non 
robust and we flag them in the figures in the paper. 

3.2. Reliability of the fit 

We test the robustness of our size estimation by ap- 
plying the same fitting procedure to a set of simulated 
galaxies. We generate 1,000 galaxies with randomly in- 
put magnitude (21 < zgso^m < 24), effective radius 
(0.1" < Re,in < 1-2"), a Sersic index following a gaus- 
sian distribution {jj,, a) = (4, 2) (with the constraint 
riin > 0.1, to prevent from negative values), random 
position angle and axis ratios following a gaussian dis- 
tribution (/x,ct) = (0.65,0.1). The magnitude and axis 
ratio ranges are representative of our sample. The range 
in Sersic index and effective radius is chosen accord- 
ing to the local distribution of ETGs (iCaon et all 119931: 
iBlanton et all [20051: IShankar et al.|[2010[ r 

We then convolve the simulated galaxy with the PSF 
image and add Poissonian noise. The simulated galaxy 
is eventually placed in a stamp extracted from the real 
HST/ ACS zs5o image, randomly chosen between ten po- 
sitions devoided of sources, thus taking into account the 
background noise and all possible systematics inherent 
to the image. 

In Figure [U we compare the estimated and input pa- 
rameters (magnitude: (zsso.out — zs5o,in), effective ra- 
dius: SRe = [Redout - Re m) I Re, in-, Sersic iudcx: 5n = 
{nout — 'nin)/nin) versus the input effective radius Re,in 
(left panel) and the measured effective radius Re,out 
(right panel). In the same figure, we bin the values on the 
X-axis: the mean value and the scatter for each bin are 



shown. For the right panel, we also show in red an his- 
togram (arbitrary units) of the Re distribution for our 
real ETGs. The green dashed lines delimit the possi- 
ble SRe values due to the range of the simulation input 
values Re,in- For example, because of the definition of 
SRe and of 0.1" < Re.in < 1.2", a galaxy with Re. out 
will necessary have a corresponding value of SRe within 
[-Re,o«t/l-2" — l,i?e,o«t/0.1" — 1]. 

When looking at 5Re as a function of Re,in (left panel), 
we observe that our method recovers the effective ra- 
dius with no significant bias, except for large galaxies 
{Re,in ^ 1"); whcrc it slightly underestimates (by ^ 5%) 
the radius, because a significant part of the light is lost in 
the background noise. When looking at SRe as a function 
of i?e,o«t (left panel), we again observe no significant bias 
except for the galaxies with Re,out > 1.2", which have 
their size overestimated. This is a direct consequence 
of the chosen range for i?e,m ([0.1", 1.2"]): as the green 
dashed line illustrates, all our simulated galaxies with 
Re,out > 1-2" can only have their size overestimated. 
Our real ETGs never have values of derived Re so high, 
as shown in the red histogram. These correlations prop- 
agate to magnitudes and Sersic indexes. 

In the range of our data (0" < Re. out < 1-2"); sizes 
and Sersic indexes are recovered with systematics smaller 
than 8% and the magnitudes with systematics smaller 
than 0.08 mag, and also with a relatively small scatter. 
Hence our estimates of magnitude. Re and n are well 
recovered, in the range covered by our observations. Us- 
ing the maximum scatter for binned simulated data, we 
assign an error of 20% to our measured Re and n. 

3.3. de Vaucouleurs vs Sersic profile 

To reduce uncertainties in the fit, we tried reduc- 
ing the numbe r of fr ee parameters by using a fixed 
Ide Vaucouleur j ()1948D profile (n = 4). We estimated 
how such a fit would be reliable as compared to a Sersic 
profile. We perform our fits again (real ETGs and simu- 
lations) with the same method, but this time with a de 
Vaucouleurs profile. We then compare the results with 
those obtained with a Sersic profile (n free) in Figure O 
The left panel represents ^^gso^out ~ zs5o,out and the right 
panel log(r"^^(/re^o„t) as a function of the output Sersic 
index Uo-ut- Real ETGs are represented with large sym- 
bols (Lynx cluster: red dots, Lynx group: blue triangles, 
CDF-S: green stars) and simulations with black dots. 
ETGs with non robust structural parameter estimates 
are plotted as empty symbols. Real ETGs measurements 
and simulations show the same trend: measuring the size 
by assuming a de Vaucouleurs profile introduces a signifi- 
cant bias, which depends on the Sersic index of the ETG. 
Those results are in qualit ative a greement with those of 
[DLOnofrio ct al. {20M) and lTavIo r et al. (2010). We wiU 
use size and n estimates obtained with a Sersic profile 
hereafter. The tables in Appendix A present sizes and 
surface brightnesses derived with both Sersic and de Vau- 
couleurs profiles. We present in Figure [TU] in Appendix B 
the Sersic index distributions for our sample. We remark 
that the presence of few ETGs with small Sersic indexes 
is not unexpected, as they have been visually selected 
through their morphology. 

4. KORMENDY RELATION (KR) 
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Figure 1. Robustness of the fitting method: comparison of the estimated and input parameters (magnitude: (zsso.out — -2850, in) i effective 
radius: SRe = {Re,out — Re., in) / Re, ini Sersic index: 5n — {f^out )/nin) versus the input effective radius Re,in {left panel) and the 

measured effective radius Re,out {right panel). In the right part of each panel, we bin the values with x-axis value bins: the plots show 
the mean value and the scatter for each bin. The red histogram in the right panel represents in arbitrary units the Re distribution for 
our real ETGs. The green dashed lines delimit the possible SRe values due to the range of the simulation input values of Re, in- because 
of the definition of SRe and of 0.1" < Re,in < 1.2", a galaxy with Re. out will necessary have a corresponding value of SRe within 
[i?e, out/1-2" — 1, -Re, out /O.l" — 1]. For 0" < Re, out < 1.2", the size range spanned by our sample (see the text for more details), there are 
no significant systematics in the Re and n parameters. Using the maximum scatter for binned simulated data, we assign an error of 20% 
to our measured Re and n. 
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Figure 2. Comparison of the estimated 2350 magnitude {left panel) and the estimated Re {right panel) when using a dc Vaucouleurs 
profile (n=4) or a Sersic profile (n free): real ETGs are represented with large symbols (Lynx cluster: red dots. Lynx group: blue triangles, 
CDF-S: green stars) and simulations with black dots. ETGs with non robust structural parameter estimates arc plotted as empty symbols. 
For both magnitude and Re, we observe a systematic bias depending on the Sersic index n. 



A powerful tool to investigate the ETG evolution and 
constra in the underlying process e s is the Fundamenta l 
Plane (jPiorgovski fc Davis! Il987t iDressler et all I1987D . 
which is a scaling relation between the effective radius 
Re, the mean surface brightness (/i)e and the central ve- 
locity dispersion (Tq • As obtaining velocity dispersions of 
ETG at z ^ 1 is observationally expensive, many stud- 
ies have focused on the projection of the Fundamental 
Plane along the velocity dispersion axis, k nown as the 
Kormendy Relation (KR) ()Kormendvlll977[) : 

{^)e = a + 13 X log Re, (3) 

where Re is in kpc. The value of a depends on the pho- 
tometric band and on the redshift. Th e slope P has been 
found to be constant out to z = 0.64 l|La Barbera et al.) 



[2OOI . 

We convert our zs5o magnitudes in the i?-band rest- 
frame, Bzo, in order to derive the _B-band rest-frame 
surface brightness {ijP)e- We use the index zq to refer 
to the rest-frame and Zobs to refer to the observed frame. 
To esti mate Bz„, we u se a method similar to the one 
used in lMei et al]()2009[) . We use CB07 models (choosing 
BC03/MA05 models changes by less than 0.1 mag) 
and consider a set of galaxies with a redshift of formation 
1.8 < Zform ^ 7, a solar metallicity and an exponentially 
declining SFH with 0.1 < SFH r (Gyr) < 1. We then hn- 
early fit the relation between the colors {Bzg — zs^q,z^^^) 
and («775,Zob, -2850,206.) (where 1775, z^^, and zs5o,z,i, are 
the apparent magnitudes in the 1775 and zg^o bands for 
galaxies observed at z = Zobs)- Once this relation is es- 
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tablished, we can estimate the total i3-band rest-frame 
magnitude Bz^ from the fuU measured apparent mag- 
nitude in the ^775 and zsso bands (pubhshed in Rll). 
Eventuahy, we transform this magnitude into mean sur- 
face brightness by averaging half of the total flux on the 
surface within Re and correct for the cosmological dim- 



Zobs 



Bz„ + 2.5 X log{2TiRl) - 10 x log{l + Zobs)- (4) 



ming (1 

Taking into account the different steps in estimating 
(/i^}e, we assign an uncertainty of 0.4 mag for our (/^^)e 
estimate. 

The KR we obtain is plotted in Figure [3] The upper 
panels show our KR for the three environments: Lynx 
cluster ETGs (left panel a), red dots), Lynx group ETGs 
(middle panel h), blue triangles) and CDF-S field ETGs 
(right panel cj, green stars). ETGs with non robust 
structural parameter estimates are plotted as empty sym- 
bols. For each environment, the colored area represents 
the 1(7 dispersion of the best linear fit to our data, done 
through a classical chi-square error statistic minimiza- 
tion. The red dotted line represents a line of constant 
absolute magnitude = —20.2 mag, corresponding to 
our cut in selection at ^850 = 24 mag. The black solid and 
dashed lines represent the local KR: they represent the 
best linear fit and its la dispersion to the data measured 
in the B-band bv lJorgensen et all ()1995D for 31 ETGs in 
the Com a cluster (converted to the AB magnitude sys- 
tem). As lJorgensen et ahl ()1995[ ) sizes are estimated with 
a de Vaucouleurs profile and we have demonstrated that 
this changes the s ize estimate (Figu r e [21), we use an ap- 
proach similar to ILa Barbera et all (j2003[) and exclude 
from this sample the three largest galaxies (log(i?e/kpc) 
> 1) for which the size difference between a de Vau- 
couleurs and a Sersic profile is likely to be significant. In 
order to show that this choice of local relation does not 
affect our conclusions on the KR, we display in Figure 
111! in Appendix B a figure similar to Figure [31 but with 
our sizes estimated with a de Vaucouleu rs profile instead 
of a S ersic profile and with including in I Jorgensen et ahl 
(|1995D local relation the three largest galaxies. 

4.1. KR: dependence on the environment 

As shown in previous works at z > 1, the 
KR is in place at z ~ 1. 3 in the field (e.g. , 
di Serego Alighieri et all l2005t iLonghetti et~all 120071: 



Cimatti et al.ll2008l : IDamianov et all 20091: Saracco et al l 
20091) and in cl usters (e.g.. IHolden et al.l 120051 



Rettura et al.l [20Tol ). We find that, though the range 
in size is similar for the three environments, the distri- 
bution of cluster ETGs seems to be more concentrated 
towards smaller sizes. We will come back to this point 
in §??. We then plot the data for our whole sample in 
the lower left panel d) and the Icr dispersion around the 
best linear fit relations in the lower middle panel e). The 
KRs in the three environments are in agreement, and we 
do not observe any d ependence of the KR, on the envi- 
ronment at z ^ 1.3. iRettura et al.l (|2010D studied the 
KR in the field and in a cluster at z ^ 1.2 and found no 
dependence on the environment. Our work confirms this 
study and extends its results to the group environment. 

4.2. KR: comparison with the local relation 



When comparing with the local KR, we observe that 
our relation is shifted towards brighter luminosities and 
the slope is steeper. This change in slope may be linked 
to the magnitude cut due to the depth of our zgso image 
(see the line showing the depth of our ACS image in 
Figure [3|), or be a real steepening of the KR. 

Stellar population models predict that the luminosity 
evolution depends on galaxy age and its SFH. In Figure[31 
we show these dependences between z ~ 1.26 and z = 0. 
We plot the luminosity evolution function of age 
at the given redshift for several exponentially declining 
SFHs with characteristic time SFH r ranging from 0.1 to 

1 Gyr. The range in SFH r encompasses the likely values 
for our ETGs: for all models (BC03/MA05/CB07), our 
estimated maximum S FH r is below 1 Gyr f or 90% of our 
sample (Rll; see also IRettura et al.l[20TTl ). If evolving 
passively down to z = 0, a 3 Gyr old ETC at z — 1.26 
will be 1.5-2.5 mag less luminous in B^^ whereas a 1 Gyr 
old ETC at z = 1.26 will be 3-3.5 mag less luminous in 
Bz„ . Older ETGs evolve less in B^^ . 

In Figure [5l we code our galaxy ages (as derived from 
Rll, see this paper for details) in gray levels, for the three 
models (BC03/MA05/CB07). Larger ETGs tend to be 
older. To better visualize this trend, we bin our data in 
three size bins (log(i?e/kpc) < 0, < log(i?e/kpc) < 0.5 
and log(i?e/kpc) > 0.5). For each size bin, we overplot 
as magenta squares the mean and the Icr dispersion for 
log(i?e/kpc) and {^^)e values. We report, in orange at 
the bottom of the figure, the mean and standard devia- 
tion of the estimated ages for each bin. Thus, we observe 
an age gradient in our KR, larger ETGs being on average 
older, because they are on average more massive (see also 
Rll; such an age grad ient is observed in t he local Uni- 
verse, see for instance iShankar et al.l [20101 ). Under the 
assumption of only passive evolution in luminosity and 
according to stellar population models, this age gradient 
should lead to a steepening of the slope of the KR with 
increasing redshift, in qualitative agreement with what 
we observe. 

4.3. KR: comparison with literature data at z 1-2 

On the panel of Figure [3l we overplot as squares with 
black outlines the KRs published for ETGs at Zsp e c ^ 1 - 

2 . The sample from Idi Serego Ahghieri et al.f (|2005f ) 
(in magenta) is composed of 16 field ETGs (from the 
K20 survey, selected according to their spectra) at 0.9 < 



■'Spec 



< 1.2 (we removed the two ETGs at z, 



0.67) 



The size is estimated by using Gim2d (jSimard et al.l 
I2002D and fitting a Sersic profile on HST/ ACS zsso- 
band and VXr/FORS-l z-band images. The sample 
from iCimatti et al.l (120081 ) (in light blue) is composed 
of 13 passive galaxies (6 in the field and 7 in a cluster- 
like structure) selected from the GMASS project, mainly 
ETGs, with 1.4 < Zspec < 2. The size is estimated us- 
ing Galfit by fitting Sersic profiles to HST /ACS zsso- 
band images. For those two samples, the radius is the 
circularized effective radius. We converted the surface 
brig htness to the AB magnitude system for the sample 
of idi Serego Ahghieri et al.l (|2005l ). 

From this comparison, we can observe two facts. 
Firstly, we observe that our KR is broadly con- 
sistent with those tw o stud ies. The KR from 
Idi Serego Alighieri et al.l (120051 ). observed at lower red- 
shifts, is slightly shifted towards fainter luminosities and 
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Figure 3. Upper panels a), h), c): Kormendy relation for our ETGs at 2 ~ 1.3 in rest-frame _B-band. Lynx cluster ETGs {a), red 
dots), Lynx group ETGs (b), blue triangles) and CDF-S ETGs (c), green stars). ETGs with non robust structural parameter estimates are 
plotted as empty symbols. For each environment, the colored area represent the la dispersion of the best linear fit to o ur data. Typical 
uncertainties are represented by the cross in the lower left corner. The black solid and dashed lines represent the local KR llJorgensen et aO 
1199511 and the la dispersion for the Coma cluster, respectively. The red dotted line represents a line of constant absolute magnitude 
= —20.2 mag, corresponding to our cut in selection at zg^o = 24 mag. Lower panels d), e), f): Kormendy relation for the three 
environments simultaneously. In (d) we show our data and in (e) the areas corre sponding to the la scatte r of a linear fit. In (/): same as 
panel e), but with data from the literature; dSA05 corresponds to the sample of Idi Serego Alighieri et al.l I I2005I1 and COS corresponds to 
the sample of lCimatti et al] lj2008l ) (see text). Our Kormendy relation does not depend on the environment and is in qualitative agreement 
with passive evolution when comparing with 2 ~ and z ~ 1-2 literature data. 
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the KR from iCimatti et aLl (|2008[ ). observed at higher 
redshifts, is lying on the higher luminosity side of our 
KR. Thus, putting together those three KRs, we see a 
shift of the KR towards bright luminosities with increas- 
ing redshift, qualitatively consistent with passive lumi- 
nosity evolution. 

Secondly, looki ng at the range in size , we notice 
that the sample of Idi Serego Alighieri et all (|2005[ ) lacks 
galaxies smaller than 1 kpc, even if comparable to our 
sample in a Ks-hand limit magnitude, and is com- 
parabl e to our samp l e for large galaxies. The sam- 
ple of ICimatti et al.l (|2008( ) lacks galaxies larger than 



~3 kpc and is comparable to our sample for small 
sizes. We remark that the observed lack of small/large 
galaxies in those two samples is not a selection ef- 
fect due to the depth of the images, which would pro- 
duce a cut along a line parallel to the red dotted line 
(|di Serego Alighieri et al.f (j2005| ) limiting magnitude is 
about = -20.1 mag and ICimatti et al.l ([2008i ) data 
are deep enough to detect a galaxy at = —20.2 mag). 
Our sample ranges a larger interval in size that both 
other samples. 




log (Re/kpc) 

Figure 5. Kormendy relation at z ~ 1.3: dependence on ago. Our KR is plotted with ages estimated from tliree models 
(BC03/MA05/CB07) whicii are coded in gray levels. ETGs with non robust structural parameter estimates are plotted without black 
outlines. The solid and dashed black lines and the red dotted lines are the same as in Figure |3] We bin our data in three size bins (vertical 
dashed magenta lines): for each bin, we overplot as magenta squares the mean and the la dispersion for log(iJe/kpc) and {^i^)e values. 
We also report, in orange at the bottom of the figure, the mean and standard deviation of the estimated ages for each bin. On average, 
large ETGs are older than smaller ones. 



5. MASS-SIZE RELATION (MSR) 

In Figure |6l we plot our MSR, derived using the 
three stellar population models (BC03/MA05/CB07), 
and splitting our sample by environment: we display 
from upper to lower panels, Lynx cluster ETGs (red 
dots), Lynx group ETGs (blue triangles), GDF-S ETGs 
(green stars) and all environments simultaneously. The 
solid and dashed lines represent the local MSR scaled 
to a Salpeter IMF, and its la relatio n, resp ectively. The 
local MSR estabhshed bv lShen et all ()2003l ) with sizes es- 
timated in z band for SDSS galaxies selected according 
to their Scrsic index (n > 2.5) is in cyan. The local MSR 
established bv I Valentinuzzi et al.l (j2010al ) with sizes esti- 
mated in V band for WINGS cluster galaxies that were 
morphologically selected to be ETGs is in magenta. The 
difference in the rest-frame used to estimate sizes would 
shift the local MSRs towards larger sizes (around 10%- 
15% accordi ng to iBcrnardi ct al. 200 31). We compare our 
sample with I Valentinuzzi et al.l (|2010al ). because both se- 
lect ETGs from a morphologicall y classificat i on. A ll our 
results do not change when using IShen et al.l (l2003l) local 
MSR, which is widely used in the literature. 

5.1. MSR: dependence on the environment 

We plot in Figure [3 for the three environments and 
the three models, the normalized distributions of the size 
ratio Re/ Re,Vaien., which represents the ratio between 
the size of our ETGs and t he one predicted by the local 
MSR of lValentinuzzi et al.l (l2010al ) at similar masses. For 
each histogram, we overplot with a black solid line the 
best-fit gaussian to the distribution, obtained through a 
non-linear least-squares fit. As found in previous studies 
at z ~ 1-2, our sample presents a significant number of 
ETGs having small radii compared to the local ones of 
similar mass. However, the precise number of such ETGs 
and the value of size ratios depend on the model (and the 
local MSR used as a reference) . We display in Table [1] 
the mean and standard deviation corresponding to the 
gaussian fit. For qualitative comparison, we also display 
in this Tab le the corre s pondi ng values when comparing 
our sizes to lShen et all (|2003D local MSR. 

Despite the dependence on the model (and on the lo- 
cal MSR), there is a general trend : most of the clus- 



ter and group ETGs lie below local MSRs. We observe 
in Figure [7] that cluster and group ETG size ratios are 
mostly below 1 with a narrow distribution that peaks 
around 0.6-0.8, whereas field ETG size ratios have a 
more widespread distribution, peaking around 0.7-1.1. 
The values in Table [T] confirm this point, i.e. that at a 
given mass, ETGs in denser environments tend to have 
smaller sizes at z ^ 1.3 than in the local Universe. From 
a Kolmogorov-Smirnov (Kuiper) statistical test, the clus- 
ter and field samples for ETGs with masses M < IO^^Mq 
do not have the same size ratio distributions, at 85% 
(90%) and 90% (95%) using MA05 and CB07 models, 
respectively. Using BC03 stellar population models, on 
the other hand, the null hypothesis (the cluster and field 
samples are taken from the same statistical distribution) 
cannot be rejected (rejected at only 40% and 60% con- 
fidence for a Kolmogorov-Smirnov and Kuiper test, re- 
spectively), in this mass range, and it is rejected at 86% 
and 90%, respectively, on the entire mass range. We un- 
derline that the Kuiper test is more sensitive to the shape 
of the distribution than the Kolmogorov-Smirnov test. 

We remark that our size ratios in clusters are in agree- 
ment with previous estimates in high reds hift clusters 
(|Rettura et all 120101 : IStrazzullo et all l2010f ). We know 
that ETGs with emission lines have hi gher size ratios 
'e^ lToft etaLl [2007t iZirm et al.l [2007t I Williams et al 



201C). Even when we remove them (|Salimbeni et al 



20091 Holden, Nakata, private communication) our dis- 



tributions remain similar as in Figure [71 as it can be seen 
in Figure [121 in Appendix B. 

The field ETG population is approximately equally di- 
vided with ETGs above and below the local MSRs, ex- 
cept when using the BC03 models, in which case a ma- 
jority lie below the local MSRs. This can be explained by 
the underestimate of the TP-AGB phas e by BG03 mod- 
els. A s explained in previous works (e.g. iMaraston et al.l 
[20061 [Conroy & Gunn 2010, Rll), this underestimate 
leads to an artificial increase of the estimated age and 
mass of galaxies with ages ~l-2 Gyr (see Figure 5 & 6 
of Rll). Such an ETG will be estimated with an older 
age and a greater mass: it will be shifted towards the 
more massive area in Figure [6l and thus will more likely 
lie below the local MSRs. This effect of BC03 models is 
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Figure 6. Mass-size relation at 2 ~ 1.3 derived witii tiie three stellar population models (BC03/MA05/CB07), split by environments: from 
upper to lower panels, Lynx cluster ETGs (red dots), Lynx group ETGs (blue triangles), CDF-S ETGs (green stars) and all environments 
simultaneously. ETGs with non robust st ructura l parameter estimates are plotted as empt y symbols. Magenta (resp. cyan) solid and 
dashed lines represent the local relation of I Valent inuzzi et al. (2010a) (resp. [5hen et aT]|2003l ) and its \a dispersion. Typical uncertainties 
are represented by the cross in the right left corner. When quantifying size evolution, one has to be careful to the model and local MSR 
used (see Table [TJ. Most of cluster and group ETGs lie below local MSRs, whereas field ETGs are in agreement with local MSRs (except 
with BC03 models, for which more field ETG lie below the local MSRs, see text for discussion). 



Table 1 

Size Ratios as a Function of Stellar Population Model, Local MSR, and Environment 



Environment 


BC03 


Re 1 Fte^V alen 

MA05 


CB07 


BC03 


Re/Re, Shen. 
MA05 


CB07 


CLE-W 


0.7 ± 0.3 


0.6 ± 0.2 


0.8 ± 0.2 


0.5 ± 0.2 


0.5 ± 0.1 


0.7 ± 0.2 


GRl-3 


0.6 ± 0.2 


0.8 ± 0.2 


0.8 ± 0.2 


0.6 ± 0.2 


0.7 ± 0.2 


0.8 ± 0.2 


CDF-S 


0.8 ± 0.5 


1.1 ± 0.7 


1.0 ± 0.6 


0.7 ± 0.3 


1.1 ± 0.6 


1.0 ± 0.5 



All 



0.7 ± 0.4 



0.7 ± 0.2 



0.8 ± 0.4 



0.6 ± 0.3 



0.6 ± 0.2 



0.8 ± 0.3 
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Figure 7. Size ratio Re/ Re,v alen normalized distributions derived with the three stellar population models (BC03/MA05/CB07), split 
by environments: Lynx cluster ETGs {upper panels, red tilted lines), Lynx group ETGs {middle panels, blue tilted lines) and CDF-S 
ETGs {lower panels, green horizontal lines). Re/Re Valen represents the ratio between the size of our ETGs and the one predicted by the 
local MSR of I Valentinuzzi et al.l Il2010al ). The black solid line represents the best-fit gaussian to the distributions. The black dashed line 
represents the locus of Re /Re Valen = 1- For MA05/CB07 models, field ETG population is in agreement with the local MSR, whereas 
cluster /group ETG population lie below the local MSR. 
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more obvious for our CDF-S sample, because this sample 
contains more ETGs with ages ~l-2 Gyr (see Rll). 

In Rll, we underlined that our CDF-S sample might 
be biased against low- mass/passive ETGs, because their 
low luminosity and the lack of emission line prevent to 
derive reliable spectroscopic redshift. We have tested 
that our conclusions for the field sample do not depend 
on this potential bias. We have selected a new sample of 
GOODS /CDF-S ET Gs, using photometric redshifts from 
iSantini et al.l (|2009( ) and we obtain a MSR that again 
shows a distribution similar to the local. The selection 
criteria are described in Appendix B and the results are 
in Figure [T51 

If on the other hand, we are missing massive galaxies, 
this would not change the overall mass-size distribution, 
which clearly shows to be similar to the local one at all 
masses and will not be changed significantly by rare mas- 
sive galaxies. 

5.2. Size ratio versus redshift of formation/stellar mass 

We now check the dependence of the size ratio on the 
redshift of formation. We plot in Figure [5] the size ra- 
tio Re/ Re,Vaien. as a fuuctiou of the redshift of forma- 
tion Zform- For indication, we also mark with thick 
black outline the ETGs known to have emission lines 
and with a thick orange outline the ETGs known to 
be passive (jSalimbeni et al.l [20091 Holden, Nakata, pri- 
vate commun ication). As already observed in the lit- 
eratu re (e.g. iValentinuzzi et afl l2010at iWilliams et al.l 
I2O1O0 . when looking at our resuhs with MA05/CB07 
models, we observe that passive/quiescent ETGs tend 
to have small size ratios, whereas line-emitting/ star- 
forming ETGs tend to have larger size ratios. 

If we consider the plots with MA05/CB07 models (for 
which the redshift of formation is more robust because of 
the better modeling of the TP- AGB stellar phase) , most 
of cluster and group ETGS lie below the local MSRs, 
as already shown in Section 5.1. We also observe that 
galaxies with small size ratios do not have a preferred 
redshift of formation and, on the other hand, there is 
a deficit of ETGs with z/orm ^ 3 and a high size ratio 

{Re/ Re.Valen. ^ 2). 

Figure [9] shows the dependence of the size ratio on 
the stellar mass. We do not observe any clear depen- 
dence of the size ratio on the stellar mass. In particu- 
lar, for MA05/CB07 models and all environments, ETGs 
with small size ratio span the whole range in mass of 
our sample. Moreover, the high-mass end of our sample 
(M ^ 2 X 1O^^M0, see also Figure |6l) is in agreement 
with the local MSR for all environments. 

5.3. MSR: comparison with literature data at z 1-2 

Our results are consistent with published works of the 
MSR at the same redshift and with comparable mass 
r ange, i.e. IQi" < M / M^ < 3 x 10". 

IValentinuzzi et all (|2010b[ ). whose sample covers 
masses slightly higher to those of this work up to a red- 
shift z ~ 0.7, have found that the median MSR in galaxy 
clusters have only mildly evolved between z ^ 0.7 to the 
present (size ratio ~ 0.5-0. 8 when c o mpar ing ETC pop- 
ulations). Working with a iKroupal (j2001| ) IMF and on 
a mass-selected sample [Mxroupa > 4 x lO^^Mo), and 
defining a superdense galaxy as a galaxy with E50 — 



a.bxAlKr 



— -— > 3 X 10^ A/0 kpc ^, these authors find that 

41% of their z ^ 0.7 cluster sample is made of super- 
dense galaxies, against 1 7% for z ~ cluster sample 
(jValentinuzzi et ani2010aD . Using the very same criteria, 
we find that ~ 45% (BC03) or - 35% (MA05/CB07) of 
our cluster sample is made of superdense galaxies. The 
mass cut reduces our cluster sample to ^ 15 galaxies, 
lowering the statistics: assuming Poissonian errors on 
our number of galaxies leads to an uncertainty of ^ 20% 
on our estimated percentages. A correct estimate should 
also take into account the different galaxy selection and 
si ze estimate met h ods. 

iNewman et al.l (|2010[ ) study a sample of field 
spheroidals with masses 6 x 10^" < M/Mq < 3 x 10" 
in the redshift range 1.1 < Zgpec < 1-6. They find that 
galaxies less massive than IO^^Mq lie on the local MSR, 
and those more massive than 10^^ Mq have to grow of 
around twice in size. Those results are consistent with 
o ur founding in t he fie ld. 

iSaracco et al.l (j2011| ) study field and cluster (among 
which Lynx members) ETGs lying at 0.9 < z^pec < 2. 
They find that compact ETGs formed over a wide range 
of redshift (2 < Zform < 10) and that normal ETGs have 
formed at z < 3. Our results are in agreement, since 
we find that our cluster ETGs (more compact) formed 
over a wide range of redshifts and field ETGs (that have 
younger ages) have a larger dispersion in sizes. In our 
sample these different distributions seems to be linked to 
the environment. They also find compact ETGs through- 
out all their probed mass range (5 x 10^" < AI/AIq < 
5 X 10"). 

Ivan der Wei et all ()2008( ) studied a sample of morpho- 
logically selected ETGs in field and cluster environments 
at 0.8 < Zspec < 1-2 and, using dynamical masses, found 
that ETGs have on average increase their size by a fac- 
tor of 2 between z ~ and z ~ 1. This result is at ap- 
parent discrepancy with ours for field ETGs. However, 
these two works probe different ranges in galaxy mass: 
our sam ple spans masses of 10^° < M/Mq < 3 x 10^^, 
whereas Ivan der Wei et aLl (j2008l ) sample includes more 
massive, 8 x 10^° < M/Mq < 10^^ (when dynamical 
mass es are converted to Salpet er stellar masses). In addi- 
tion, Ivan der Wei et all (|2008D sizes are estimated using 
de Vaucouleurs profile, which complicates a possible com- 
parison for the few ETGs in common with similar mass 
see ji l3.3l and Figure (H). The fact that Ivan der Wei et al.l 
2008f ) do not find an environmental dependence on the 
evolution of the MSR and we do, probably reflects the 
different range in mass probed by the two works. This 
points to a very interesting situation since it suggests 
that the evolution of the MSR is mass dependent. 

Recent works have found very com pact galax- 
ies i n the field at Zsv ec ^ 1.5 (e.g. IPaddi et al.l 

2OOI Ivan Dokkum et all [20081 : iDamianov et all 120091 : 
van Dokkum et al.l 120101 ). The mass range covered in 
these works (M > 2 x IO^A'/q, when converted to 
Salpeter stellar masses) hardly overlaps our mass range, 
and we thus cannot compare conclusions directly because 
we are sampling different ranges in mass and in redshift. 

Our results are driven by galaxies with masses M < 
2 X 10"Mq. Our galaxies with masses M ~ IO^Mq fol- 
low the same trends that the entire sample: field galaxies 
lie on the local MSR relation, cluster galaxies show an 
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Figure 8. Size ratio versus redshift of formation, derived with the three stellar population models (BC03/MA05/CB07): Lynx cluster 
ETGs are red disks, Lynx group ETGs are blue triangles and CDF-S ETGs are green stars. ETGs with non robust structural parameter 
estimates are plotted as empty symbols. We also mark with thick black outline the ETGs known to have emission lines and with thick 
orange outline the ETGs known to be passive. The black dashed line represents the locus of Re/ Re.Valen = 1- Fof MA05/CB07 models 
and cluster/group environments, ETGs with small size ratios have equally distributed z/orm and ETGs with high size ratios are younger 
ETGs {zform J; 3). For the field, our sample do not show old galaxies (see also Rll) probably because it does not cover large areas. 
However, again, field galaxies cover a larger distribution in sizes. 
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Figure 9. Size ratio versus stellar mass, derived with the three stellar population models (BC03/MA05/CB07): symbols are as in Figure 
\E\ We do not observe any clear dependence of the size ration on the stellar mass. 



average MSR shifted to sizes 30-50% smaller. Our galax- 
ies with masses M > 2 x 10^"'^ M0 are a few, and show 
a large dispersion in size: they lie on the local MSR in- 
dependently of the environment (see Figure [6l), but their 
small number does not permit us to draw conclusions on 
their behaviour. 

6. CONCLUSIONS 

In this work, we have studied a sample of 76 ETGs 
spanning a wide range of environments (cluster, group, 
and field) at z ^ 1.3, combining multi- wavelength ob- 
servations of the Lynx supercluster, with data on the 
GOODS/CDF-S field. We estimated the size of our 
ETGs by fitting a Sersic profile to the HST/ ACS zgso 
images, which probe the rest-frame B-band. Combining 
those sizes with stellar masses and stellar population ages 
derived in Rll, we are able to study two crucial struc- 
tural relations, the Kormendy relation and the mass-size 
relation, in three different environments at z ^ 1.3. 

We obtain the following results: 

1. The Kormendy relation, in place at z 1.3, does 
not depend on the environment. We thus confirm 
the result of iRettura et al.l ()2010D and extend it to 
the group environment. Our results are in agree- 



ment with results in the c l uster and field samples of 
Idi Serego Alighieri et al.l (|2005() and lCimatti et al.] 
(I2008n . 



2. Concerning the mass-size relation, for all stellar 
population m odels (BC03/MA05/CB07) and lo- 
cal rel ations CShen et al.l 120031 : IValentinuzzi et aTl 
l2010af ). ETGs are on average more compact in 
denser environments. When comparing the MSR 
at high redshift with the one in the local universe, 
the uncertainty on the mass coming from the model 
used to estimate it and the choice of the local 
MSR can significantly influence the conclusion on 
the importance of the size evolution. When us- 
ing MA05/CB07 models, we find that the major- 
ity of cluster and group ETGs are below the local 
relations, whereas field ETGs follow a MSR simi- 
lar to the local one. From a Kolmogorov-Smirnov 
(Kuiper) statistical test, the cluster and field sam- 
ples for galaxies with masses M < 10^^ Mq do 
not follow the same size ratio distribution, at 85% 
(90%) and 90% (95%) using MA05 and CB07 mod- 
els, respectively. When using BC03 models, the 
two distributions do not differ. 
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3. When using MA05/CB07 models, we find that 
compact ETGs do not have a preferred redshift of 
formation. Those results are in close agreement 
with those of Saracco et al.l ()2011h . who studied a 
sample of 62 ETGs with 0.9 < Zgpec < 2. As con- 
cluded by these authors, the lack of dependence 
of the compactness on the redshift of formation is 
not consistent with models that predict compact 
galaxies to have formed at earlier times, when the 
Universe was more dense. 

When we compare the MSR of cluster and group ETGs 
vs field ETGs, we find that, at similar masses, cluster 
and group ETGs arc more compact than field ETGs. 
On average this does not depend on cluster galaxy age. 
This result is in contrast with what has been found so 
far for field galaxies a t z 1 at higher masses (e.g. 
Ivan der Wei et ani2008| ). and it might be due to the dif- 
ferent range in masses that we are probing. If this was 
confirmed by larger samples, it would mean that envi- 
ronmental effects are visible in the evolution of the MSR 
for ETGs with M <2x 10"Mq. 

Our results are mainly driven by galaxies with masses 
M < 2 X lO"A/0. Our galaxies with masses M ^ 
IO^^Mq follow the same trends that the entire sample. 
Our galaxies with masses M > 2 x 10^ "'^Mq are a few, 
but they lie on the local MSR independently of the en- 
vironment (see Figure [6]) ; however their small number 
does not permit us to draw conclusio ns on their be- 
haviour. As concluded by other authors (iNewman et al.l 
I2010t iCassata et al.ll201lUSaracco et al.ll2011[ ). the verv 
compact galaxies at z 2 should have gone a dramatic 
evolution in size to reproduce our results at z ^ 1.3. 
This growth between z ^ 2 and z ~ 1 see ms to be some- 
how d ifferent in cluster and field galaxies. ICassata et al.l 
(|2011| ) have shown that, in the field, ETGs enlarge their 
size and increase their stellar mass by a factor of 5 be- 
tween z = 2 and z~l. Atz^l.3, field galaxies 
are already on the local MSR, while cluster galaxi es still 
have compact sizes on average (see also IStrazzuU o et al.l 
120101 for a similar result at z = 1.4), indicating that 
their size distribution stil l needs to be enlarged (see also 
IValentinuzzi et al.ll2010b[ ). 

Since in the local Uni verse, the ETG M SR does not 
depend on environment (jMaltbv et al.|[20lol) . our results 
imply that an evolution in the MSR of cluster and group 
ETG size is required to explain current observations, 
while field ETGs show a MSR that is compatible with 
the local one. The evolution of the MSR in dense envi- 
ronments might reflect either an evolution in size of the 
pristine population or the transformation of ETG pro- 
genitors that arc not classified as ETG at z ~ 1.3 or the 
accretion of a new population of larger ETGs. 

In the first case, minor dry merger events could have 
enlarged the size of the ETG population. In the second 
case, compact ETGs might have not had much evolution, 
but a new population of larger ETGs could have been 
formed by not-ETG pr ogenitors or accreted in de nse en- 
vironments at z < 1 (jValentinuzzi et al.l l2010"B l. This 
new population might not have been observed at z 1.3 
because its progenitors are not ETGs at that time. These 
galaxies might be disk galaxies that have evolved from 
a large bulg e spiral population or galaxy merger s (e.g., 
lPostman~et al.. ,2005^ ,Mei ct al., i2006ak .Poggianti et al.l 



[200llValentinuzzi et al.ll2010bl : [MeTet al.ll201lD . For in- 
stance, according to semi-analytic models, the ETG pop- 
ulation at z ~ 1 in dense environments contain less than 
^ 70 % of the stellar mas s which ends up in ETGs at 
z - (jKaviraj et all 120091 ). 

On the other hand, these results pose some challenges 
to current state-of-the-art galaxy evolution models that 
predict a nearly m ass-independent size for ETGs (e.g. 
iShankar et al .1120 111 ) and, we have checked, nearly inde- 
pendent of environment. More detailed theoretical work 
is required to fully understand all the processes at work 
that can affect galaxy sizes. This is clearly beyond the 
scope of the present work and will be the subject of fu- 
ture efforts. 
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APPENDIX 

APPENDIX A: LYNX AND CDF-S ETC STRUCTURAL PARAMETERS 
Lynx (resp. CDF-S) ETG structural parameters are presented in Tables[2]and[3] (resp. Table S} 



ETGs at z ~ 1.3. IV. Scaling relations in different environments 



Table 2 

Lynx cluster ETG structural parameters and surface brightness. 
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18.2 
















4 


0.21 


1.79 


18.8 


1486 


08 


48 


31.72 


+44 


54 


42.95 


6.8 


0.14 


1.13 


17.7 
















4 


0.09 


0.78 


16.9 


1794 


08 


48 


32.78 


+44 


54 


07.22 


3.1 


0.17 


1.43 


17.8 
















4 


0.21 


1.72 


18.2 


1922 


08 


48 


32.99 


+44 


53 


46.69 


2.5 


0.14 


1.20 


16.6 
















4 


0.18 


1.49 


17.0 


1525 


08 


48 


33.01 


+44 


55 


11.92 


8.1 


0.43 


3.61 


19.5 
















4 


0.20 


1.65 


17.8 


1962 


08 


48 


33.04 


+44 


53 


39.75 


5.4 


0.14 


1.18 


17.5 
















4 


0.05 


0.44 


15.4 


2094 


08 


48 


34.08 


+44 


53 


32.32 


2.7 


0.28 


2.31 


18.8 
















4 


0.39 


3.26 


19.5 


2343 


08 


48 


35.98 


+44 


53 


36.12 


3.4 


1.13 


9.47 


20.6 
















4 


1.44 


12.00 


21.1 


2195 


08 


48 


36.17 


+44 


54 


17.30 


6.6 


0.96 


7.99 


20.5 
















4 


0.38 


3.19 


18.5 


2571 


08 


48 


37.08 


+44 


53 


34.05 


4.0 


0.25 


2.06 


18.3 
















4 


0.25 


2.08 


18.3 



Note. — Re denotes the circularized effective radius. Uncertainties on Re 
are of 20%, uncertainties on n are of 20% and uncertainties on {fi^)e are of 0.4 
magnitudes. We did not report parameters considered as non robust. 
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Table 3 

Lynx group ETG structural parameters and surface brightness. 



ID R.A. DEC. n Be Re {lJ.^)e 

(J2000) {J2000) (arcsce) (kpc) (AB mag.arcsec-^) 



Lynx Group 1 ((z) = 1.262) 



518 


08 49 03.52 


+44 53 21.62 


4.7 


0.24 


1.99 


19.1 








4 


0.21 


1.72 


18.8 


1339 


08 49 08.32 


+44 53 48.32 


4.5 


0.54 


4.48 


19.6 








4 


0.46 


3.87 


19.2 


1024 


08 49 09.00 


+44 52 44.08 


3.8 


0.08 


0.71 


16.4 








4 


0.09 


0.72 


16.4 


825 


08 49 11.24 


+44 51 29.19 


4.5 


0.22 


1.87 


17.3 








4 


0.20 


1.70 


17.1 


1249 


08 49 12.27 


+44 52 13.05 










1085 


08 49 13.69 


+44 51 18.82 


2.7 


0.23 


1.94 


18.0 








4 


0.31 


2.62 


18.6 


Lynx Group 2 {(z) = 1.260) 


1636 


08 49 00.92 


+44 58 49.15 


3.3 


0.31 


2.58 


18.8 








4 


0.37 


3.11 


19.2 


1383 


08 49 03.99 


+44 57 23.37 


3.2 


0.19 


1.60 


18.2 








4 


0.23 


1.91 


18.5 


2000 


08 49 07.15 


+44 57 52.04 


5.6 


1.06 


8.83 


21.1 








4 


0.58 


4.83 


19.8 


Lynx Group 3 {{z) = 1.263) 


137 


08 48 53.26 


+44 44 22.39 










542 


08 48 55.14 


+44 44 58.83 










1135 


08 48 56.28 


+44 46 45.62 


5.9 


1.13 


9.45 


21.5 








4 


0.62 


5.15 


20.2 


889 


08 48 56.63 


+44 45 39.90 


1.5 


0.15 


1.25 


18.5 








4 


0.25 


2.08 


19.6 


1431 


08 48 57.31 


+44 47 08.01 


4.7 


0.23 


1.95 


18.9 








4 


0.20 


1.67 


18.6 


1064 


08 48 57.79 


+44 45 57.51 


1.9 


0.13 


1.05 


17.8 








4 


0.19 


1.55 


18.7 


1136 


08 48 57.96 


+44 46 04.53 


5.9 


0.48 


4.00 


20.0 








4 


0.30 


2.53 


19.0 


1775 


08 49 01.62 


+44 46 28.23 


6.6 


1.00 


8.32 


21.8 








4 


0.43 


3.62 


20.0 


1731 


08 49 04.43 


+44 45 08.65 


5.7 


0.67 


5.60 


20.1 








4 


0.41 


3.40 


19.1 



Note. — Re denotes the circularized effective radius. Uncertainties on Re 
are of 20%, uncertainties on n are of 20% and uncertainties on (/i^)e are of 0.4 
magnitudes. We did not report parameters considered as non robust. 



APPENDIX B: FURTHER TESTS AND EXPLANATIONS 

This appendix shows some of the tests explained in the text. 

Figure [TOl show the Sersic index distribution for our sample (see also lMei et al.l[20Tlh . 

Figure [TT] com pares our sizes derived with a de Vaucouleurs profile to the local KR. derived from the 31 ETGs 
measurements of iJorgensen et ahl l)1995f ). We show that our results from Section ?? do not change. 

Figure [12] is similar to the figure [71 but we removed from the sample the ETGs known to have emission lines and as 
a consequence of it, larger sizes (see Section 5.1). It shows that even if we do remove these galaxies our results do not 
change. 

Figure [T3| shows how the MSR in the field would change if we add to our CDF-S spectroscopic sample ETGs that 
are selected using photometric redshifts. To perform this test, we add to our spectroscopic sample a sample of ETGs 
selected to have with 1 < Zphot < 1-5 and no reliable Zspec (jSantini et al.|[2009t ). in the same magnitude and color range 
that our spectroscopic sample We thus include any possible low- mass/passive ETGs which may be absent from our 
CDF-S sample. We remark that this test sample might be contaminated by outliers. For the new ETGs included in 
the sample, we derive masses and sizes with the same procedure used in this paper. Figure [T51 shows the MSR for this 
test sample (upper panels) and the size ratio Re/ Re,Vaien normalized distributions (lower panels). The distribution 
of the size ratios Re/ Re,Vaien is similar to the one we find in Figure [T] Our results do not change: the field MSR has 
still a distribution similar to the local relation. In the main body of the paper, we leave results obtained using only 
field spectroscopically confirmed members, since we expect that the photometric redshift selected sample in the field 
might present a much higher contamination than in clusters and groups. 



ETGs at z ~ 1.3. IV. Scaling relations in different environments 



Table 4 

CDF-S ETG structural parameters and surface brightness. 



ID " 


Zspec 




R.A. 


DEC. 


n 


Re 


Re 








(J2000) 


(J2000) 




(arcsec) 


(kpc) 


(AB mag. arcsec ^) 


3680 


1.119 


03 


32 20.28 


-27 52 33.01 


4.1 


0.07 


0.61 


15.1 












4 


0.07 


0.61 


15.1 


10069 


1.119 


03 


32 19.36 


-27 47 16.24 


2.9 


0.13 


1.04 


17.7 












4 


0.16 


1.27 


18.1 


7237 


1.123 


03 


32 45.14 


-27 49 39.95 


6.5 


0.19 


1.52 


16.8 












4 


0.12 


0.99 


15.9 


3000 


1.125 


03 


32 23.60 


-27 53 06.35 


6.5 


0.31 


2.57 


18.4 












4 


0.18 


1.49 


17.2 


7567 


1.158 


03 


32 23.28 


-27 49 26.07 


5.0 


0.15 


1.25 


17.3 












4 


0.13 


1.05 


16.9 


10717 


1.173 


03 


32 30.83 


-27 46 48.56 


3.2 


0.17 


1.39 


17.8 












4 


0.21 


1.71 


18.3 


14747 


1.178 


03 


32 39.17 


-27 43 29.02 


3.6 


0.34 


2.85 


19.7 












4 


0.39 


3.24 


19.9 


9066 


1.188 


03 


32 33.06 


-27 48 07.54 


4.3 


0.34 


2.79 


19.4 












4 


0.31 


2.56 


19.2 


4176 


1.189 


03 


32 24.98 


-27 52 08.63 


4.3 


0.22 


1.85 


18.8 












4 


0.21 


1.74 


18.6 


14953 


1.215 


03 


32 25.98 


-27 43 18.93 


5.3 


0.20 


1.67 


17.8 












4 


0.16 


1.30 


17.3 


11062 


1.220 


03 


32 46.34 


-27 46 32.00 


4.4 


0.50 


4.18 


20.7 












4 


0.45 


3.72 


20.4 


15093 


1.222 


03 


32 35.63 


-27 43 10.14 


6.3 


0.42 


3.45 


18.6 












4 


0.24 


1.96 


17.4 


12264 


1.222 


03 


32 26.29 


-27 45 36.19 


1.3 


0.07 


0.55 


14.6 












4 


0.09 


0.74 


15.3 


12000 


1.222 


03 


32 26.26 


-27 45 50.71 


5.9 


0.31 


2.59 


18.6 












4 


0.21 


1.73 


17.7 


9702 


1.223 


03 


32 35.79 


-27 47 34.76 


5.0 


0.26 


2.18 


19.2 












4 


0.21 


1.74 


18.8 


4981 


1.253 


03 


32 44.26 


-27 51 26.75 


6.9 


0.30 


2.52 


18.5 












4 


0.17 


1.39 


17.2 


288 


1.264 


03 


32 25.40 


-27 56 09.88 


5.1 


0.06 


0.51 


16.0 












4 


0.06 


0.47 


15.8 


10650 


1.277 


03 


32 08.37 


-27 46 51.21 


8.2 


0.22 


1.82 


17.5 












4 


0.11 


0.92 


16.1 


6791 


1.297 


03 


32 50.19 


-27 50 01.04 










9369 


1.297 


03 


32 16.02 


-27 47 50.00 


7.7 


0.62 


5.21 


20.1 












4 


0.24 


2.00 


18.0 


10231 


1.317 


03 


32 39.63 


-27 47 09.12 


5.4 


0.92 


7.74 


20.8 












4 


0.58 


4.86 


19.8 


17506 


1.328 


03 


32 20.08 


-27 41 06.75 


6.1 


0.42 


3.50 


18.6 












4 


0.25 


2.07 


17.5 


1857 


1.345 


03 


32 38.37 


-27 54 08.83 


5.2 


0.42 


3.56 


19.4 












4 


0.31 


2.61 


18.8 


969 


1.346 


03 


32 35.99 


-27 55 09.49 


7.0 


0.55 


4.61 


20.2 












4 


0.26 


2.15 


18.5 


10041 


1.356 


03 


32 25.04 


-27 47 18.20 


0.6 


0.13 


1.12 


17.6 












4 


0.29 


2.41 


19.2 


12505 


1.374 


03 


32 06.81 


-27 45 24.35 


3.8 


0.14 


1.15 


16.9 












4 


0.14 


1.18 


17.0 


8938 


1.382 


03 


32 33.98 


-27 48 14.69 


1.4 


0.08 


0.70 


16.6 












4 


0.12 


0.99 


17.4 



Note. — " ID refers to the GOODS-MUSIC v2 catalogue of ISantini eFall l(200gf l. Re 
denotes the circularized effective radius. Uncertainties on Re are of 20%, uncertainties on n 
are of 20% and uncertainties on {fJ.^)e are of 0.4 magnitudes. We did not report parameters 
considered as non robust. 
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Figure 10. Sersic index distributions for our sample: red histogram is for cluster ETGs, blue histogram is for group ETGs and green 
histogram is for field ETGs. The dashed line represents a Sersic index of 4, corresponding to a de Vaucouleurs profile. We also report the 
mean and standard deviation for the three distributions. 




log (Re/kpc) 

Figure 11. Same as Figure [S] but witli our sizes estimated with a do Vaucouleurs profile. Moreover, we include here in the local KR of 
IJorgensen et al] I I1995I) the three largest galaxies that we excluded in Figure [S] Wc observe that none of our conclusions of Section ?? is 
affected if we use de Vaucouleurs sizes. 
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Figure 12. Size ratio Re/ Re.Valen normalized distributions derived with the three stellar population models (BC03/MA05/CB07), split 
by environments: same figure as Figure [7] but we here remove ETGs known to have emission lines (7 ETGs in clusters, 4 ETGs in groups 
and 14 ETGs in the field). Our results do not change. 
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Figure 13. Results when considering a CDF-S sample including ETGs with photometric redshift. Upper panels: MSR, symbols 
are as in Figure |6] Lower panels: size ratio Re/Re^Valen normalized distributions derived with the three stellar population models 
(BC03/MA05/CB07), same as Figure[7l This sample shows a MSR similar to the local, as the spectroscopically selected sample. This test 
shows that our spectroscopically selected sample is not biased with respect to a photometric redshift selection. 



